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We develop a semiclassical kinetic theory for electron spin relaxation in semiconductors. Our 
approach accounts for elastic as well as inelastic scattering and treats Elliott- Yafet and motional- 
narrowing processes, such as D'yakonov- Perel' and variable g- factor processes, on an equal footing. 
Focusing on small spin polarizations and small momentum transfer scattering, we derive, starting 
from the full quantum kinetic equations, a Fokker-Planck equation for the electron spin polarization. 
We then construct, using a rigorous multiple time scale approach, a Bloch equation for the macro- 
scopic (fc-averaged) spin polarization on the long time scale, where the spin polarization decays. 
Spin-conserving energy relaxation and diffusion, which occur on a fast time scale, after the initial 
spin polarization has been injected, are incorporated and shown to give rise to a weight function 
which defines the energy averages required for the calculation of the spin relaxation tensor in the 
Bloch equation. Our approach provides an intuitive way to conceptualize the dynamics of the spin 
polarization in terms of a "test" spin polarization which scatters off "field" particles (electrons, im- 
purities, phonons). To illustrate our approach, we calculate for a quantum well the spin lifetime at 
temperatures and densities where electron-electron and electron-impurity scattering dominate. The 
spin lifetimes are non-monotonic functions of temperature and density. Our results show that at 
electron densities and temperatures, where the cross-over from the non-degenerate to the degenerate 
regime occurs, spin lifetimes are particularly long. 

PACS numbers: 72.25.Rb, 72.25.Dc, 72.25.-b 



I. INTRODUCTION 



The spin degree of freedom of an electron provides an additional variable that potentially can be used to add new 
functionality to electronic, optoelectronic and magnetoelectronic devices or to even build radically new devices entirely 
based on the coherence of electron spin states. This has led to the newly emerging field of spintronics. 0,0 A subclass 
of spintronics device concepts relies on the capability to inject, control, and detect electron spin polarizations in non- 
magnetic semiconductors. j3|,t4] The spin polarization, which would enable the device operation, is a non-equilibrium 
state and its characterization, e.g., in terms of lifetimes and transport coefhcients, has to be given within a kinetic 
theory. 

Of particular interest is the lifetime of the non-equilibrium spin polarization in non-magnetic n-type III-V semicon- 
ductors. Important spin relaxation processes for itinerant electrons in this class of materials include the Elliott- Yafet 
(EY) process ||5|,Q, that leads to spin- flip scattering and, in materials without inversion symmetry, the D'yakonov- 
Perel' (DP) process in which spin states precess because of spin off-diagonal Hamiltonian matrix elements resulting 
from a combination of spin-orbit coupling and inversion asymmetry. An external magnetic field, in many cases re- 
quired to control and manipulate the electron spin, can also influence the electron spin dynamics. It quenches the DP 
process thereby tending to extend the spin lifetimes as a function of magnetic field, and it opens a spin relaxation 
channel due to the fc-dependence of the electron g-factor, which forces the spin of electrons in different quantum states 
to precess around an external magnetic field with different rates. ITOi] For brevity we will refer to this mechanism 
as a variable g-factor (VG) process. 

Spin dynamics in semiconductors has been extensively studied in magneto-optics pdL IT^ IT^ Il4j using various 
spin-sensitive emission, transmission, and reflection spectroscopies. These spectroscopies are now re adily adaptable 
to spatially and time resolved measurements jlllillHlIllialMIIlllllMlMlMlMllllMIM 
which, together with the emergence of spintronics concepts, inspired new theoretical investigations in bulk 
mHalllEi and dimension-reduced semiconductors. MEpjlMEililElEllillil&S, 49, 50, 51j,|5l|5S 

The theoretical investigations are based on the early work 0,0,0 augmented by modern band structure theory for 
bulk and dimension-reduced semiconductors. The EY spin relaxation rates are usually calculated using the Golden 
rule for spin-flip scattering, whereas the spin-flip rates due to motional-narrowing (DP and VG) processes are, at 
least conceptually, obtained from a semiclassical Boltzmann-type equation for the non-equilibrium spin polarization, 
although not always is the Boltzmann equation explicitly solved. Instead, a common procedure is to adapt the 
expression for the spin relaxation rate originally derived by D'yakonov and Perel' to the scattering processes under 
consideration. 

The D'yakonov-Perel' expression for the spin relaxation rate, which results from the solution of the Boltzmann 
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equation within the elastic approximation, treats all scattering processes on-shell, even inelastic scattering events, 
e.g., due to electron-electron or electron-phonon scattering. The obtained spin relaxation rates are therefore on-shell 
rates, which have to be averaged over energy with an appropriate weight function before a comparison with experiments 
can be attempted. Usually, the difference of the distribution functions for spin-up and spin-down electrons divided 
by the total number of electrons contributing to the initial spin polarization is used as a weight function. T^, A(L 
The ad-hoc energy averaging, which is necessary because of the incomplete treatment of inelastic scattering processes, 
can lead to substantial deviations from the spin relaxation rates obtained, e.g., from a full numerical solution of the 
Boltzmann equation. 35] 

We develop in this paper a systematic kinetic theory for electron spin relaxation, applicable to spin-flip (EY) and 
motional- narrowing (DP and VG) spin relaxation processes in bulk and quantum wells, which avoids the ad-hoc energy 
averaging and gives a clear physical picture of the time evolution of the optically or electrically injected non-equilibrium 
spin polarization. We derive, starting from the full quantum kinetic equations for the electron Green functions, a 
semiclassical Fokker-Planck equation for the time evolution of the non-equilibrium spin polarization, valid for small 
spin polarizations and for small momentum transfer scattering, and employ a multiple time scale perturbation approach 
to separate the fast spin-conserving from the slow spin non-conserving time evolution. As a result, we obtain on the 
time scale of spin relaxation a Bloch equation for the macroscopic (/c-averaged) non-equilibrium spin polarization, 
which is the quantity measured in, e.g., time-resolved Faraday and Kerr rotation experiments. [T^ ll^ ITsL IT^ I29ll30l | 
The weight function defining the energy averages needed, e.g., for the calculation of the spin relaxation tensor and the 
spin relaxation rates, turns out to be directly related to the quasi-stationary spin polarization, which is the terminating 
state of the initial, fast spin-conserving time evolution of the injected spin polarization. Our approach treats spin-flip 
(EY) and motional-narrowing (DP and VG) processes on an equal footing. Due to the different angle dependences, 
a Matthiessen-type rule holds, however, for isotropic semiconductors, where the total spin relaxation tensor is simply 
the sum of the individual spin relaxation tensors. The diagonal elements of the spin relaxation tensor, the spin 
relaxation rates, are either given in terms of an energy averaged spin-flip rate (EY process) or an energy average of 
a generalized relaxation time, which accounts for both on- and off-shell scattering events, multiplied by a precession 
rate (DP and VG processes). 

In the next section we introduce a generic model for electrons in n-type III-V semiconductors applicable to bulk 
and quantum well situations. In Sec. Ill we give a complete description of our semiclassical kinetic theory for the 
electron spin dynamics. As far as the formal development is concerned, we treat EY, DP and VG processes on an 
equal footing and also allow for quenching effects due to orbital motion of electrons in an external magnetic field. 
In Sec. IV we apply our approach to the particular situation of DP spin relaxation in an idealized quantum well at 
temperatures and densities for which electron-electron and electron-impurity scattering dominate. Our main findings 
are summarized in Sec. V. Technical details concerning the calculation of the quantum well collision integrals due to 
electron-electron and electron-impurity scattering are relegated to an Appendix. 

II. MODEL HAMILTONIAN 

We consider conduction band (CB) electrons in III-V semiconductors, e.g., GaAs, in the presence of an applied 
magnetic field. The model used here applies to both bulk and quantum well situations. Within an envelope function 
approach ,54. .55j . which treats the two states at the conduction band minimum explicitly and includes a large set of 
states perturbatively, the effective mass Hamiltonian for the CB electrons can be cast into the form 

where K = k — {e/hc)A{r) and A(r) is the vector potential. The "spin basis" for the CB electrons used to define the 
model is a = 4- and a = —, where a ~ + (a = — ) denotes a state which is mostly spin- up (spin-down) with a 
small admixture of spin-down (spin-up) . 

The first term denotes the dispersion of the Kramers degenerate conduction band which, depending on the sophisti- 
cation of the envelope function approach, could contain nonparabolicity effects. For quantum wells k and e(fc) denote 
the in-plane momentum and the in-plane dispersion of the conduction subband under consideration. The second term 
comprises the Larmor precession due to the external magnetic field, with HCIl ~ 1^b9* B the Larmor energy vector. 
Here /i^ and g* denote the Bohr magneton and the electron g-factor. The third term describes spin off-diagonal 
Hamiltonian matrix elements arising from the coupling to higher lying states. The most important of which are the 
splitting of the conduction band due to inversion asymmetry (lA) and the term which leads to a fc-dependent electron 
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g-factor. For bulk semiconductors, the two contributions are given by 

hftiAiK) = 2Jo«/A(i?) , (2) 
hng{K) = 2aiK'^B + 2a5{K,B-K} + 2a6T{K,B) , (3) 

respectively. The definition of the vectors kia{K) and f{K,B) and of the parameters Sq and Ui can be found in 
Refs. [53.l55j | and {.., ..} denotes an anticommutator. The expressions for conduction subbands in a quantum well are 
obtained by averaging the bulk expressions 121 and Q over the subband envelope function. 

In addition to bulk inversion asymmetry, dimension-reduced semiconductors can have additional sources of asym- 
metry due to interfaces which share no common atom ;3l5 or due to layer design (structural inversion asymmetry (39j). 
Both mechanisms can be cast into spin off-diagonal Hamiltonian matrix elements and can be therefore treated in the 
same way as the spin off-diagonal terms due to bulk inversion asymmetry. 

For a complete description, a collision term arising from electron-impurity, elect ron-phonon and electron-electron 
scattering, 

B^c ^ei ^" Hf3p -\- , 

is added to the effective mass Hamiltonian. The electron-impurity term reads 

- E E ^'^"."'(^' ^')cLc£.„. , (4) 

kk' 

with a scattering matrix element given by 

3 

The Bloch states for the conduction band are not pure spin states, because of spin-orbit coupling. The scattering 
matrix element contains therefore an overlap factor 

/oa'(fc,fc') = R,fc|/7„,,fc,) , (6) 

which is of order unity for a = a' (spin conserving scattering) and is small, but not zero otherwise (spin non-conserving 
scattering). The electron-phonon collision term would have the same structure as Eq. Q but with phonon creation 
and annihilation operators appearing in the matrix element Maa' {k,k'). The electron-electron scattering contribution 
has the form 

TJpp = — M„, (/ci , fco, fcs, ^4)ct cl cr c,- , (7) 

2 /-^ aia2Q3Q4V U ^; 1/ j,^^^ j,^^^ ksOs k^at ' V I 

kiOLi 

where the scattering matrix element, 

MQ,jQ2a3a4(fci, fc2, fca, /C4) = V {ki " fc4)/Qi (fci , fc4)/a2a3 (/C2 , fcs)^^^ ^ _^ , (8) 

contains two overlap factors. The functions U{k) and V{k) denote, respectively, the potential of a single impurity 
(neutral or ionized) and the Coulomb potential between two conduction electrons. 

The model Hamiltonian is characterized by e(fc), r2/^(fc), ^lg{k), and Ia.a'{k,k'). These quantities need to be 
obtained by an electronic structure calculation. The formal structure of the kinetic theory described in the next 
section is independent of the particular form of these quantities. 



III. SEMICLASSICAL KINETIC THEORY 



In this section we give a systematic derivation of the Fokker-Planck equation governing the electron spin relaxation in 
the limit of small spin polarizations. The derivation is independent of dimensionality, applying to bulk semiconductors 
and semiconductor heterostructures, and treats motional-narrowing (DP and VG) and spin-flip (EY) spin relaxation 
processes on an equal footing. To obtain a Fokker-Planck equation, we restrict ourselves to the Born approximation, 
but collective effects giving rise to dynamical screening of the Coulomb interaction can be approximately incorporated 
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at the level of a quantum analog to the Lenard-Balescu equation. |5y] Besides its intuitive interpretation in terms of 
a small "test" spin polarization scattering off a bath of "field" particles (impurities, electrons and phonons), causing 
dynamical friction, diffusion, and eventually relaxation for the "test" spin polarization, the Fokker-Planck equation 
is the starting point for a multiple time scale analysis which results in the derivation of a Bloch equation for the 
macroscopic (fc-averaged) spin polarization. Its decay is usually characterized by the diagonal elements of a spin 
relaxation tensor, which are quadratures of either a spin-flip rate (EY process) or a generalized relaxation time 
multiplied by a precession rate (DP and VG processes). 

Since the derivation is quite lengthy and to some extent rather formal we first give a short outline of the main steps. 
We start from the full quantum kinetic equations for the Keldysh Green functions, [s^, IH, Each component of 

the Keldysh Green function is a 2 x 2 matrix in electron spin space. In the first step we derive, within the semiclassical 
approximation, a kinetic equation for the density matrix. This accounts to treating momentum scattering processes 
as instantaneous on the time scale of spin relaxation, which is usually the case. Calculating the self-energies which 
appear in the semiclassical kinetic equation in the Born approximation, linearizing with respect to spin polarization, 
and expanding the self-energies up to second order in the momentum transfer (diffusion approximation) finally yields 
a Fokker-Planck equation for the spin polarization, which we then analyze in terms of multiple time scale perturbation 
theory. 



A. Kinetic equations 

For a spatially homogeneous system (we assume a constant magnetic field B) , the information about spin relaxation 
is contained in the electronic density matrix, which, due to the spin degree of freedom, is a 2 x 2 matrix in spin space, 

N^,c.Ak,t) =< ici^^Cf:J{t) > , (9) 

but diagonal in fc-space. Here, the operators evolve in time with the full Hamiltonian, including the time-dependent 
perturbation, which could be, e.g., a circularly polarized light pulse applied at time t — tQ. To perform the averaging 
in Eq. (jSJ denoted by < [...] >, we consider the system to be in thermodynamical equilibrium for t < tg, take the 
limit to ~oo and evaluate the expectation value in Eq. Q with respect to the equilibrium density matrix. I6C 

To derive a kinetic equation for the density matrix it is convenient to start from Keldysh Green functions. 
For a constant magnetic field, the vector potential is a function of f. It is therefore necessary to initially work with 
kinetic equations in real space. In this subsection we set h = 1. Introducing a numerical index 1 that stands for 
fiaiti and 2 for r2Ct2t2, we write in the notation of Ref. |59| 

^G^2 = * ( g|+ g|' ) ' ^^^^ 

Note that each component of the Keldysh Green function is a 2 x 2 matrix in spin space. Introducing further a 
self- energy 




y++ y+- 
^12 ^12 

^12 ^12 



(11) 



we set up two matrix Dyson equations, one where the time differentiation is with respect to ti and one where it is 
with respect to 12'- 

dt,Gi2 = -iTA2 - «(eG)i2 - ir,(EG)i2 , (12) 
dt,Gi2 = iT,6i2 + «(Ge)i2 + i{Gt)i2T, , (13) 

with (5i2 — 5{ti — t2)5{12) and the energy matrix ei2 = 5{ti — t2)i{l2) = 5{ti — t2)'?o^(12)e(— iV^^ — (e/c)A(ri)), 
where we neglect nonparabolicities in the dispersion. We adopt the convention that numerical indices written as a 
subscript contain the time variable, whereas numerical indices written as an argument do not. Matrix multiplication 
with respect to the Keldysh indices is implied and internal variables are summed (integrated) over; Tz is a Pauli matrix 
and To is the unit matrix in Keldysh space. 
Subtracting Eq. from Eq. gives 



[L,G]i2 = f,(EG)i2-(GE)i2f, , 



(14) 
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where [.., ..] denotes the commutator. To compactify the notation, we introduced a differential operator 

Li3 = fo(5i3i(3) = foSi3{idt3 - e(3)) . (15) 

It is understood that in the second term of the commutator, the operator L32 acts to the left with the temporal 
differential operator dt^ replaced by its adjoint —dt^- 

Equation (|14|l contains two time variables. To obtain a kinetic equation for the electronic density matrix, which 
depends only on a single time variable, it is necessary to perform the equal time limit. This is most conveniently done 
in the (mixed) Wigner representation, where the equal time limit reduces to an integration. Separating the self-energy 
into a singular and a regular part [53 

Sf^ = AP^{12;h)5p,6ih - h) + S?! , (16) 

introducing relative and center variables, r — ri — r2, R — (fi + ^2)/2, t — ti — t2, and T — {ti + 12)/2, and defining 
a Fourier transformation with respect to the relative variables. 



= exp 1 



A{R,T,k,oj)^ dr dre"^^-'''^A{R,T,r,T) , (17) 

J —00 J 

together with a gradient operator [g^I 

' ' ■95-9:^a# + v4.v|-V|.Vf] , (18) 
the equal time limit of the ++ component of Eq. H14() can be written as 

DiR, T, k) = F{R, T, k) + C{R, T, k) , (19) 

with a driving term on the Ihs, 

/OO J 
{g^^L{R, T, k, uj)G++{R, T, k, Lu) 
-00 27r 

- g^^G++{R,T,k,oj)L{R,T,k,iu)) , (20) 
and a rhs, which contains a molecular field term arising from the singular part of the self-energy, 

/OC 7 
-00 27r 

_ ^GAg++(^^j.^^^^)^++(^^y^^)] ^ (21) 

and a colhsion term due to the regular part 

-00 27r 

+ g^^±+-{R,T,k,Lj)G-+{R,T,k,uj) 

- g^^G+-iR,T,k,uj)f:-+{R,T,k,uj)) . (22) 

The semiclassical approximation amounts to the assumption that the Green functions and self-energies vary slowly 
on the macroscopic scales, T and R, respectively. It is therefore sufficient to keep in Eq. p9|) only the leading order 
terms in a gradient expansion. The leading order of the rhs of Eq. H19() is the zeroth order, i.e., Eqs. (|21|l and H22() 
with g^^ 1. The Ihs of Eq. p9|) . however, has to be determined to first order, because the zeroth order vanishes. 
Using 

L{R,T,Lu,k)=uj + e{K) , (23) 

we explicitly obtain 

/OO T 
^ [dTG++ {R, T, fc, c^) - e{K) ■ V^- G++ (i?, T, k, u;) 

+ V£ e(i?) • G++ (i?, T, fc, u;) . (24) 
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To ensure gauge invariance of the kinetic equation we follow Ref. |6l| and consider the generalized momentum 
K ^ k — {e/c)A{R) as an independent variable instead of the momentum k. Using the two identities [Glj . 

a{K) = a{K) ■ V + a(K) x (V^ x K) , 
d-bx{Vj^xc) ~ {a -V c) -h — {b -V c) ■ a , 

Eq. becomes 

D{R,T,K)^(dT + ^Vge{K)-{BxVj^)+Vge{K)-V^)i j ^G++iR,T,K,Lu) . (25) 

We assume weak interactions and replace the full Green functions by the noninteracting Green functions (quasi- 
particle Ansatz), 

GP^iR, T, K, Lo) = Gl\k ^ K, c^)U„(£)^^(^,^,;f) , (26) 

where the noninteracting density matrix N(){k) is replaced by the full density matrix N{R,T, K). Performing the 
cj— integrations in Eqs. H21|l . H22I) . and (|25() then yields the semiclassical kinetic equation for the electronic density 
matrix: 

(^T + ^ e{K) -BxV^ + Vji e{K) ■ V^) N{R,T, K) = 
i[N{R, T, k), A++(i?, T, k) + ±++{R, T, K, e{k))] 

+iN{R,T,k)t-+{R,T,k,e{k)) + it+-{R,T,k,t{k))[l^ N{R,T,k)] . (27) 

To obtain a closed kinetic equation for the electronic density matrix, internal Green functions, which appear in 
the self-energies, have to be of course also eliminated according to Eq. (|26|l . Details concerning the calculation of 
self-energies are given in the next section. 

For a homogeneous magnetic field, the electronic density matrix does not explicitly depend on R. The R dependence 
can be therefore neglected. For a quadratic dispersion, e{K) — K"^ /2m* (?i = 1 in this subsection), the Lorentz term 
moreover becomes 

(e/c)V^ e{k) X B • V^, (i? x ^c) ■ = -i^c ■ C , (28) 

where ^Ic — eB /m*c is the cyclotron energy vector and £ the angular momentum operator in iiT-space, and we obtain 
the kinetic equation for the electronic density matrix in a more familiar form: 

(dr-inc- tjN{T,K) = i[N{T,k),l\++{T,K) + t++{T,k,e{K))] 

+ iN{T,k)Y,-+{T,k,e{k)) 

+ zE+-(r,i?,e(/?))[l-iV(T,i?)] . (29) 

This equation is the basis for the calculation of the spin relaxation time in spatially homogeneous systems subject 
to a constant magnetic field. The first term on the rhs describes the coherent motion in a molecular field modified 
by correlation effects. If S^^, and N where scalar functions, as in ordinary transport theory, this term would 

vanish. The molecular field term is therefore a consequence of the quantum mechanical treatment of the spin degree 
of freedom. To the singular part of the self-energy contribute the spin off-diagonal terms in the Hamiltonian and the 
Hartree-Fock fields due to electron-electron scattering. Dissipation and relaxation originate from the regular part of 
the self-energy and give rise to the second and third term on the rhs. They are at least second order in the interaction. 
Formally, they correspond to the scattering-out and the scattering-in terms in a matrix-Boltzmann equation. The 
matrix structure is of course a consequence of the full quantum mechanical description of the spin. Only momentum 
scattering is treated classically. 



B. Calculation of the self-energies 

The semiclassical approach to furnish the self-energies in the matrix-Boltzmann equation H29|l . valid for magnetic 
fields, which do not restructure the electron dispersion, is to represent interaction processes in terms of diagrams, 




(e) J2) (f) 

A A Q Cl 

(g) 



FIG. 1: Diagrammatic representation of self-energies in the Born approximation for electron-electron (a)-(d), electron-impurity 
(e), and electron-phonon scattering (f). Diagram (g) denotes the self-energy due to spin off-diagonal Hamiltonian matrix 
elements. 

calculate the diagrams using standard rules [s^ to obtain i, t') and A^''(fc, i), perform the zeroth order gradient 

expansion, and then replace all internal Green functions according to Eq. H26() . This heuristic strategy leads to self- 
energies, which can be expressed in a manifestly gauge invariant form by writing the internal momentum integrations 
in terms of the generalized momentum K. The formal structure of the self-energies is then the same as without 
magnetic field. Within the semiclassical approach, the magnetic field gives therefore only rise to the Lorentz term. 

In Fig. ^we depict the self-energies in the Born approximation for electron-electron (a-d), electron- impurity (c), 
and electron-phonon (f) scattering. Diagram (g) corresponds to the self-energy due to the spin off-diagonal term in 
the Hamiltonian 

The Hartree-Fock diagrams (a) and (b) contribute to the instantaneous self-energy A++. They are second order 
in the spin polarization and therefore, for sufficiently small spin polarizations, negligible. There are two second order 
diagrams due to electron-electron scattering, the direct (c) and the exchange (d) Born diagram. Anticipating that soft 
scattering dominates, we neglect the exchange diagram (d). The direct Born diagram contributes to 5]++, and 
S ^ . It can be shown, however, that E++ is at least second order in the spin polarization and therefore negligible 
in the limit of small spin polarizations. The E and components, contributing to the collision integral, are in 
contrast linear in the spin polarization and cannot be neglected. Diagrams (e) and (f), corresponding to the Born 
approximation for electron-impurity and electron-phonon scattering, only contribute to E^"?. As in the case of electron- 
electron scattering, the H — V component can be again neglected, if the spin polarization is small enough, whereas the 
H — and — + components contribute in linear order in the spin polarization to the collision integral. Diagram (g), 
corresponding to the spin off-diagonal part of the Hamiltonian is linear in the spin polarization and contributes 
to Eventually it leads to a torque force acting on the spin polarization. 

Anticipating small spin polarizations, we neglect E++ and the Hartree-Fock contribution to A++. Writing in the 
Born approximation furthermore 1b [A^] = I'b \^\ + ^'b + \^\ fo'" second and third term on the rhs of Eq. 
(|29|1 . the semiclassical kinetic equation for the electronic density matrix reduces to 

[hdt - ihnc ■ tj N{k, i) = ^ N{kt), (hClL + hCliAik) + fiflg(fc)) -S +Ib[N] , (30) 

where we relabeled the center time T ^ t and adjusted to the notation of Eq. ^\ h \s explicitly included and the 
generalized momentum is now denoted by k. 

Equation H3U|) is a matrix Boltzmann equation similar to the semiconductor Bloch equations frequently used to 
describe optically pumped semiconductors. [60( Thus, numerical techniques used for the solution of the semiconductor 
Bloch equations can be adopted to the numerical solution of Eq. (|30|l . Calculations of this kind have been successfully 
performed for various situations. |3^ Isslli^. IslL Is^ Issj 

To avoid a numerical solution, we focus on small spin polarizations and linearize the Born collision integral with 
respect to the spin polarization. It is important to note that the equilibrium density matrix N(,q{k) = N(k,t oo) 
is not diagonal in the "spin basis" . Expanding the equilibrium density matix in terms of Pauli matrices yields 



N,g{k) = ,f{k) + -a ■ S,g{k) , (31) 



where f{k) = {l/2)TiNeq{k) = {f+{k) + /_(fc))/2 is half of the sum of the equilibrium distribution functions of the 
spin up and spin down electrons and Seq{k) is the equilibrium spin polarization. Accordingly, we also write for the 
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density matrix at arbitrary times 



Nik, t) = f{k) + Sfik, t) + -a- [Segik) + SS{k, t)] , (32) 



with 6f{k,t) and 5S{k,t) the changes induced by optical pumping or by electrical injection. We defined Neq{k) for 

t — + oo, that is, it contains the electrons created by the perturbation and both 6S{k, t) and Sf{k, t) have to vanish for 
t — > oo. 

Inserting the expansion (|32|l into the Boltzmann equation (|3U|I yields two kinetic equations, one for the charge 
component Sf(k,t) and one for the spin component 6S{k,t). The collision terms couple the two equations. If 
however only a small portion of the total number of electrons initially contributed to the spin polarization, i.e. if 
6f{k,t) <C f{k), the coupling can be ignored and it suffices to focus on the equation for 6S{k,t) alone. 

Since the total spin polarization Seq{k) + dS{k,t) is small, we linearize the Born collision integral with respect to 
both Seq{k) and 6S{k,t). Thus, 

Ib [N] = Ib [N,q] +lB[f + 6f, SS] . (33) 

If we now apply Tr(T[...] on both sides of Eq. (|Sn|) . use InWeq] = as well as CljA{k) x Seq{k) = 0, because, by 
construction, the equilibrium density matrix Ng^q{k) commutes with iJ/^i, and ignore furthermore Sf in /_b[/ + ^/: '^'5'], 
we get a closed kinetic equation for the non-equilibrium spin polarization 

(hdt - ihiic ■ -C) 5S{k, t) = (^hilL + nniAik) + nflg(fc)) x 6S{k, t) + Jb[/, 5S] (34) 

with Jb[/,<55] =Tra/B[/,<55]. 

The collision integral can be further simplified if we split the spin-flip matrix into a leading spin conserving diagonal 
part and a small off-diagonal part which describes spin-flip scattering. Since nonparabolicities are small, the diagonal 
part is approximately equal to the unit matrix and we get /(fc, fc') ~ 1 + 51{k, k'), with (5/(fc, k') <C 1. Expanding the 
collision integrals up to second order in (5/(fc, k') gives 

Js[/,55] = Jb[/,,5S^](")+ Jb[/,<5S^](i) + Jb[./,<55](2) . (35) 

The first order term J^^'' [/, 5S] potentially mixes EY, DP, and VG spin relaxation channels, but for semiconductors 
with high symmetry it does not contribute to the spin relaxation rates. 
The kinetic equation for the excess spin polarization becomes therefore 

{hdt - ihiic ■ -C) 6S{k, t) = hilL X 6S{k, t) + {hiiiA{k) + hVlg{k)^ X 6S{k, t) 

+ 4\f^5S] + fi\f,5S] . (36) 

This equation contains motional-narrowing (DP and VG) and spin-flip (EY) spin relaxation processes on an equal 

footing. The Elliott- Yafet process is simply encoded in ' [/, 6S] whereas the motional-narrowing processes result 
from the combined action of the torque forces given by the second term on the rhs and the spin conserving scattering 

processes comprising 5S]. 

Independent of the scattering process, the structure of the collision integrals in Eq. H3()|l is {y — ei, ee, and ep) 

4"' [/, 5S] = ^ [VF^^ (fc + q; q) 5S{k + q, t) - W {k; q) SS{k, t)] , (37) 
q 

[/, SS] = 2j2W''{k + q; q) g{k,k + q)X [g{k, k + q}xSS{k + q, t)] , (38) 

where, for concise notation, we introduced a spin-flip vector 

^ /Im/+_(fc,fc')\ 
9{k,k)= I Re/+_(fc,fc') I , (39) 



with (fc, fc ) the off-diagonal element of the overlap matrix ijHJi. This is a result of the Born approximation and 

the linearization with respect to the spin polarization. In general, the structure of the collision integrals depends on 
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the scattering process. Here, however, the scattering process enters only through W'^{k;q), the probabiHties for a 
transition between momentum state k ^ q and k. For electron-ionized-impurity scattering, for instance, 

^^^,0) = 2TTN,\U{q)\H{e{k - q) - e{k)) , (40) 

while for electron-electron scattering, 

W^^k;q) = 47r|T/(9)l'E{[l-/(^-'?^-/(^'+'?^]/(^') 

k' 

+ f{k-q)f{k +q)^5{e{k)^e{k)-e{k-q)-e{k +q)) , (41) 

with U{q) and V{q) statically screened Coulomb potentials. Similar expressions hold for electron-phonon scattering. 
For electron impurity scattering, which is elastic, W'^'^{k + q;q) = W^'^{k; q); moreover iy^*(fc; q) is independent of the 
equilibrium distribution of the spin-up and the spin-down electrons. In general, however, the transition probabilities 
depend on the equilibrium distribution of the electrons, and, in the case of electron-phonon scattering, also on the 
equilibrium distribution of the phonons. 



C. Diffusion approximation 

The simple form of the collision integrals H37|l and (|38|l suggests to conceptualize the dynamics of the non-equilibrium 
spin polarization in terms of spin-polarized "test" electrons, scattering off an equilibrated bath of "field" particles 
(impurities, electrons, and phonons). Usually this picture can be only applied to electron-impurity and electron- 
phonon scattering, where the scattering partners belong to different species, and not to electron-electron scattering, 
where the scattering partners belong to the same species. It is only within the linearized spin dynamics, which 
essentially treats the electrons comprising the non-equilibrium spin polarization as a separate species, that the "test- 
field-particle concept" can be applied to electron-electron scattering as well. We now take full advantage of the 
simplicity of the collision integrals and expand the collision integrals with respect to the momentum transfer q. As a 
result the integro-differential equation becomes a differential equation. 

The on-shell spin-conserving process due to elastic electron-impurity scattering yields 

4^'°'V,SS] = E E ...gfc^ ^S{k,t) , (42) 

whereas the inelastic spin- conserving processes due to electron-electron or electron-phonon scattering give rise to an 
off-shell contribution 

nn 

4"'^°«[/,<55]= E E E a,. ...dh ^^■■■■■-■.(^)'^^(^'^)' (43) 

CC,Cp 71 

where, in both cases, the moments are defined by {u — ei, ee, ep) 

C^,...,.„(fc) = ^Y.1^^-1^r.W{k;q) . (44) 
q 

The transition probability W'^{k;q) depends on the precise modeling of the elementary scattering process and also 
on the dimensionality of the system. In the Appendix we give explicit expressions for electron-electron and electron- 
impurity scattering in a quantum well. Note, for inelastic scattering the differential operators act on the moments 
Cf^ i„ (^) whereas for elastic scattering the moments are in front of the differential operators. 

Equations 1)42(1 and 1(43(1 involve partial differential operators of arbitrary order. To obtain tractable equations, the 
expansion is in many cases truncated after the second order term (diffusion approximation). As a result, scattering 
processes with small momentum transfer are treated exactly whereas scattering processes with large momentum 
transfer are treated approximately. Because the transition probability for the (unscreened) Coulomb potential diverges 
for small momentum transfer, soft Coulomb scattering events dominate, and the diffusion approximation is expected to 
describe Coulomb scattering reasonably well. A similar reasoning applies also to electron-LO-phonon scattering. The 
singular behavior of the collision integrals is stronger in three than in two dimensions. The diffusion approximation 
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is therefore somewhat better justified for bulk than for quantum well situations. |38l 1461 1481 . |62| Nevertheless, our 
numerical results suggest that even for a quantum well the diffusion approximation gives reasonable quantitative 
results for the spin relaxation time. 

Keeping therefore only the second order terms, we write 



4"' [f- SS] - SS] + /°^-°^[f, 5S] 



dkidkj 



E 



1/— cc.cp 

V{k)dS{k,t) , 



dkidkj 



SS{k,t) 

SS{k, t) 



(45) 



where the first two terms on the rhs come from elastic scattering processes and the last two terms encode inelastic 
scattering events. In Eq. (|45|l we introduced for the first and second moments {v = ei,ee,ep). 



m.Wik-q) , 



(46) 
(47) 



which have the meaning of fc-dependent dynamical friction and diffusion coefficients, respectively. Within the diffusion 
approximation the spin-conserving (Born) collision integrals are therefore represented by a Fokker-Planck differential 
operator (|45|) . Each scattering process gives rise to a particular Fokker-Planck operator, with particular dynamical 
friction and diffusion coefficients. 

In the same spirit, expanding the spin- flip collision integral Jg [/, SS] up to second order in the momentum transfer 
q, and using g{k, k) ~ 0, gives 



with a spin-flip tensor 



-Ilik)SS{kt) , 



R{k) = 4 E^«'^(^)G^^(^) . 



(48) 



(49) 



i^— ci,cc,cp IJ 



given in terms of the total diffusion coefficient and a tensor Gy (k) which describes the rate of change of the spin-flip 
vector g{k, k'): 



with 



Gy(fc) 



Of^ + 01] 



_d_ , 



k'=k 



dkjkk' 



k'=k 



(50) 



(51) 



In Fig 121 we illustrate the physical content of the diffusion approximation encoded in Eqs. (|45|l and 148|l : the small 
"test" spin polarization SS{k,t) scatters off equilibrated "field" particles, which, depending on the scattering process, 
are either electrons, phonons, or impurities. Spin conserving scattering can be elastic and inelastic, because the "field" 
particles can absorb/emit energy, the bath has "internal degrees". Spin non-conserving scattering, on the other hand, 
turns out to be elastic within the diffusion approximation. 

We now introduce scaled atomic units and measure energy in units of a scaled atomic Rydberg Rq — Rq/s and 
length in units of a scaled atomic Bohr radius qq — ^/saQ, with RqUq — fi,^/2mo and = 2y/sRQdQ, where uiq is the 
bare electron mass and s is a scale factor chosen to yield Rq — 1 meV. Symmetry- adapted coordinates are then a 



11 




FIG. 2: Graphical illustration of the collision terms within the diffusion approximation: A "test" spin polarization scatters 
off a generalized bath of equilibrated "field" particles (electrons, impurities, and/or phonons). Whereas spin-fiip scattering is 
elastic (on-shell) within our approximation, spin-conserving scattering can be on- or off-shell, i.e. the "test" spin polarization 
can lose/gain energy, because the "internal degrees" of the bath can absorb/emit energy. 



radial coordinate e = k'^ and a generalized angle variable which, for bulk semiconductors comprises two angles, 
the polar angle 9 and the azimuth angle 0, and for quantum wells is simply the polar angle cj). Before we express the 
Fokker-Planck equation in these symmetry adopted coordinates, we recall that the experimentally measured quantity 
is the macroscopic (fc-averaged) spin polarization. Normalizing the macroscopic spin polarization to TV,, the (small) 

number of initially spin-polarized electrons and writing the fc-integral in symmetry- adapted coordinates, we define a 
"macroscopic" spin polarization (per spin polarized electron) , 



1 



de / dcu J{e)6S{e,uj,t) , (52) 



{2TT)dns Jo 

with d the dimension, Us = Ng/L'^ the density of initially spin-polarized electrons and ,]{£) the energy dependent part 
of the Jacobian, which arises from the transformation to the symmetry-adapted coordinates. Note that due to the 
normalization, 5S{Q) is a unit vector in the direction of the initial spin polarization. For bulk, d='i, du — d(j)d6 smO, 
and J{£) — a/e/S, while for quantum wells d = 2, dcu = dcfi, and J{e) — 1/2. Instead of setting up the Fokker-Planck 
equation for 5S{e,U!,t) it is more convenient to directly construct the Fokker-Planck equation for 

SS'{e,u;,t)^-^^6S{e,u;,t) . (53) 

\^7T j Tig 

The differential operator describing spin conserving scattering processes in the Fokker-Planck equation for SS' reads 
in symmetry-adapted coordinates 

= ^(^)-i4)^V), (54) 

where the operator >C^(w) denotes the total angular momentum operator in momentum space. To obtain this generic 
form for both bulk and quantum wells it is essential to include J{e) into the definition of the spin polarization. 
The off-shell term I?(e) originates from inelastic scattering events, e.g., due to electron-electron or electron-phonon 
scattering. The relaxation rates l/T/(e) and 1/Td(e) denote the rate with which the "test" spin polarization loses 
energy and the rate with which the "test" spin polarization diffuses in energy space, respectively. The on-shell term, 
describing randomization of the angle variable, is given by the last term on the rhs of Eq. H54() . It is proportional to 
the total on-shell relaxation rate l/rj_(£) due to both elastic and inelastic scattering processes. 

The rates characterizing the differential operator 'D{k) are obtained from a direct calculation of the coefficients 
A'({k) and B'(^{k) and casting the resulting differential operator I?(fc) in the specific form given in Eq. (|54() . An 
explicit calculation of the symmetry-adapted form of the relaxation tensor R(fc), which is defined in terms of the 
total diffusion coefficient Bij{k), due to both elastic and inelastic scattering processes, shows moreover that it can 
be expressed in terms of the same scattering rates. Thus, the three scattering rates 1/t/(£), l/rd(e), and 1/T±{e) 

completely specify the two collision integrals Jg'^ [/, SS] and Jg'' [/, SS] . In the Appendix we give explicit expressions 
for the relaxation rates due to electron-electron and electron-impurity scattering in a quantum well. 
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The dimensionless, symmetry- adapted Fokker-Planck equation for SS', which is the basis for the calculation of the 
spin relaxation rates presented in the next subsection, can be therefore written as 

^6S'{e,cu,t) = Pie) - -^^C^uj) +tnc ■ C{iu)]6S'{e,uj,t) + X SS'{e,cu,t) 
at 4r^(e) 

+ [Q.iA{e,Lo) + ng{e,uj)]x5S'{e,uj,t)~'R.{e,Lj) 5S'{e,uj,t) . (55) 

It contains spin relaxation due to motional-narrowing (DP and VG processes) and spin-flip scattering (EY process). 
The former arises from the combined action of the off- and on-shell spin-conserving scattering events encoded in the 
differential operators I?(e) and {\/t^{uj))C'^{uj), respectively, and the torque forces due to f1/yi(e,ci;) and f2g(e,tj), 
while the latter originates from the spin-flip tensor R(e, lo). Note that if more than one scattering process is considered, 
the relaxation rates characterizing the differential operators and the spin-flip tensor are total relaxation rates, due 
to whatever elastic and inelastic scattering processes are included in the model. The orbital motion of the electrons, 
which leads to a quenching of the motional-narrowing spin relaxation processes, is given by the term i^lc ■ on 
the rhs of Eq. H55|) . In the next subsection we develop a scheme which separates the fast spin conserving scattering 
processes from the slow spin decay causing processes and directly yields the time evolution of the macroscopic spin 
polarization. 



D. Multiple time scale analysis 

The Fokker-Planck equation (|55|l determines the time evolution of the non-equilibrium spin polarization on the 
fast, spin-conserving time scale, where randomization of the angle variables (direction of the momentum) and energy 
relaxation and diffusion occurs, and on the long time scale, where spin non-conserving processes lead to the decay of 
the spin polarization. The two time scales are well separated. The fast, spin conserving stage, whose scale is given by 
the flrst term on the rhs of Eq. H55|) and therefore by the off- and on-shell relaxation times (as well as the time it takes 
to complete one cyclotron orbit), terminates in a quasi-stationary state, which then evolves on the time scales set by 
the Elliott- Yafet term, the torque forces due to the spin off-diagonal Hamiltonian matrix elements, and the external 
magnetic field (Larmor precession). Experimentally relevant is usually the time evolution on the long time scale. In 
this subsection we employ therefore a multiple time scale approach to extract from the Fokker-Planck equation (|55|) a 
Bloch equation, which controls the time evolution of the macroscopic (fc-averaged) spin polarization on the long time 
scale. To simplify the notation, we suppress the prime. It is understood that the spin polarization dS{e,uj,t) contains 
the factor J{e)/{2TT)'^ns. 

As a preparatory step we first consider that part of the Fokker-Planck equation (|55|1 . which is spin-conserving: 

^<55(e, cj, t) = [V{e) - ^-^C'{cj) + tCic ■ £{oj)]SS{e, cj, t) . (56) 

To find the stationary solution of Eq. H56|l we set the Ihs to zero, write 

5Sst{e,uj) ^p{e)^st{e,uj) , (57) 

with p(e) defined by 

V{e)p{e) = , (58) 
and obtain for the auxiliary vector 5'st(£,a;) the differential equation 

[^*(^) " r^'^'('^) + ■ £>)]^.t(e, c^) = , (59) 



with 



s d 9^ 

V*{e) = — -- + —-— , (60) 
Tf{e)de Td{e)de^ 



the adjoint operator to I?(e). The simplest solution of H59|l is a constant vector 



(61) 
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giving rise to a stationary solution 

SSstie,io)^pie)e. (62) 

The particular form of p(e) does not matter at this point. 

We now turn to the full Fokker-Planck equation (|55|l . The spin-conserving stage of the time evolution, described by 
the first term on the rhs, occurs on a very fast time scale and is usually experimentally not resolved. Hence, it is not 
necessary to explicitly keep track of it. Instead, it is sufficient to use the final state of the fast, spin-conserving time 
evolution, i.e. the stationary solution of the spin conserving part of the Fokker-Planck equation (viz: Eqs. (|56|l - l|62() '). 
as an initial state for the time evolution on the slow time scale, where spin decay occurs. We write the initial condition 
therefore as 

SSst{e,u,0)^p{e)e, (63) 

where e is now the direction of the initial spin polarization. This initial condition is general enough, because neither 
electrical nor optical spin injection produce anisotropic initial spin polarizations. Accordingly, we write for arbitrary 
times 

SSie,Lu,t)=pie)^{e,Lu,t) , (64) 
where the time-dependent auxiliary vector 5'(e,aj,t) satisfies now the time-dependent equation 

+ (nL + niA{e,Lu) + ng{e,uj)^ x'^{e,uj,t)^-R{e,Lu)'i{e,uj,t) , (65) 

with an initial condition ^{e, ui, 0) — e. Note that Eq. (|65|l is the adjoint Fokker-Planck equation. The function p{e) 
satisfies the homogeneous differential equation (|58() and is therefore defined only up to a normalization constant. From 
the initial condition for the macroscopic spin polarization, SS(0) = e, we conclude that p{e) has to be normalized 
according to J de J dujp{e) — 1 (recall that we redefined SS such that it contains the factor J{e)/ {2TT)'^ns). Thus, 
de J dujp{e) can be interpreted as the probability density for finding a spin-polarized "test" electron in the energy 
interval [e, e + de]. 

To proceed further we scale each term in Eq. H65() to its typical value. In the case of degenerate electrons the 
typical values would be the ones at the Fermi energy, whereas for non-degenerate electrons the typical values could 
be the ones at the average thermal energy. Denoting typical values by a caret, we introduce scaled quantities t' = t/t, 
t'j.- = Tf/if, Td/Td, t'^ = t^/t±, ^'ja = ^lATiA, % = ^g^'gj ^'l = ^ltl, R' = 'Rtr, and = QcTc- The 

rescaled equation for \I''(e, tj, i') becomes (suppressing the arguments of the various functions) 



d ^, , i e d i d"^ i 1 ^, i 



ov Tf t'j: ae Td oe-^ t± 4t^ tc 



+ Ln'^ X ^' + -J-n'iA X + ^n' X - ^R'#' . (66) 

TL TiA Tg ^ TR 

We identify three time scales. A fast time scale given by the spin conserving relaxation times fi {i = f,d,l.) and 
the time it takes to complete a cyclotron orbit fc, an intermediate time scale given by the time it takes to complete 
a precession around the intrinsic magnetic fields (due to the spin off-diagonal Hamiltonian matrix elements) tja 
and Tg, and a long time scale, on which Larmor precession and spin- flip scattering occur, fi and tr, respectively. 
For representative experimental set-ups, the typical time scale t, on which the spin polarization has to be tracked 
("observation time"), and the three typical intrinsic time scales obey the following ordering: jj, 4-, J-, = 0(77"^), 



, ^ 0{ff), and 4^, 4^ — 0{r]^) where we introduced a small parameter rj. Accordingly, we classify each term 



^ Oiff), and -^,4- 

in Eq. H66|) by the smallness parameter 77. Suppressing the primes, Eq. (|66|l is rewritten as 

^-[V* - -^C? + i?lc ■ + {?liA + l1„) X ^ -f 7/ X ^ - 77 . (67) 
at 7] 4tx 

Equation (|67|) is in a form where fast and slow processes can be clearly identified. The fast spin-conserving terms and 
the orbital motion enter in order 77"^, the precession around the internal magnetic fields enters in order 77*^, whereas 
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the Larmor precession and the spin- flip scattering terms appear in order ij^. Naturally, taking as much advantage as 
possible of the existence of the small parameter 77, the first thought is to expand all quantities with respect to rj and 
apply perturbation theory. The structure of Eq. 167(1 indicates however that regular perturbation theory will lead 
to non-uniformity in the long-time regime, i.e. precisely in that regime, which we are interested in. To obtain the 
correct long-time behavior of the solution of Eq. (|67|l a multiple time scale approach is required. 

In the spirit of multiple time scale perturbation theory 63], we consider therefore ^' as a function of three time 
variables t„ ~ rj^t, n — —1,0,1, which are assumed to be independent, and substitute a second order expansion of 
the form 

into Eq. (|67|l . where the time derivative is extended to dt = rf^dt_^ + dt^ + ?79ti . Equating coefficients of like powers 
of r\ yields an hierarchy of equations for the functions Up to 0{t]) they read 

^(0) ^ [p* - -^/:2^zflc (69) 

4r_i 

#(") = \D* - -^£2 ^ . ^ (^^^ ^ ^ \ ^ ^(0) ^ (70) 

4t^ V / 

m + (1) + ^^(0) ^ [p* _ 1 £2 ^ . ^-j^(2) + (f^^^ + ^\ X ^(1) 

dt-i dto dti 4:T± V 

+ X - R^^") . (71) 

For the analysis of this set of equations it is convenient to split into an angle averaged and a remaining part, 

(£, Lu, t_i, to, ti) = a(")(e, t_i, to, ti) + Sa^''^ (e, lu, t_i, to, ti) , (72) 

where a*^"-' =< ^'^"^ >t^, with an angle average defined by 

((...))^ = / du;{...) (73) 



d 







d 






dt_i 


d 




d 


9t„i 


^dfo 


d 


^(1) 


d 


^dfo 





and < (Ja*^"-* >lj~ by definition. Since the angle variables are periodic this partitioning is always possible. From 
the initial condition, ^(£,aj,0) = e, we infer the intitial conditions a'^"^(e, 0, 0, 0) = eSn.o, and (5a''"' (e, 0, 0, 0) — 0. 
Recalling that the factor J{e)/ {2Tr)'^ns is included in the definition of SS{£,uj,t), the macroscopic non-equilibrium 
spin polarization defined in Eq. H52|l can now be rewritten as 



SSit)^ I d£p(e)a(£,t) = (a(e,t))e , (74) 

where we defined an energy average 

((...)), = / depie){...) . (75) 



Note that the function p{e), which determines the terminating state of the fast, spin-conserving time evolution, enters 
here naturally as a weight function. Formally, the weight function appears in our theory because of the Ansatz H64|l . 
which enabled us to switch to the adjoint Fokker-Planck equation. The expansion of 5* implies an analogous expansion 
for the macroscopic (fc-averaged) spin polarization: 

6S{t) = 5S^(") (t_i, to, ti) + r7<5S^(i) (t_i, to, ti) + r^^SS^^^ (t_i, to, ti) . (76) 

We now calculate the leading order contribution SS'^^^t) uniformly valid for all times. As a result, we will obtain a 
Bloch equation which determines the long time behavior of the macroscopic spin polarization. 
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1. 0{rj ^) equation 

With the substitution H72I) . the 0{ri^^) equation H69fl sphts into two independent equations, one for the angle 
averaged part and one for the angle dependent part. 

d 



d 



a(0) = l?*a(0) , (77) 



Solutions of Eqs. H77|) and (|78|l compatible with the two initial conditions, a*^°) (e, 0, 0, 0) = e and Ja^*'^ (e, 0, 0, 0) — 0, 
are a(°)(£, io, ^i) = a(°)(io,ti), with a(°)(0,0) e, and (5a(°) (e, tj, to, ti) = 0. Using Eqs. GH) and (EHI), we 
find therefore 

5S^''\t^M)^a^°\hM) , (79) 

that is, due to our choice of the initial condition, the macroscopic zeroth order spin polarization is independent of the 
fast spin conserving time scale t_i and solely evolves on the long time scales tg and ti. 



2. 0{rf) equation 

To determine the time evolution of the macroscopic spin polarization on the long time scales and ti, we study 
the 0{rf') equation. Substituting H72|) into Eq. H7UI) . the latter splits into two independent equations: 

^ -a(i) = , (80) 



dt-i dto 

^ -da'-^^ = [V* - -^C? + i?ic ■ C]da'^^^ + (diA + (ig) x a^") , (81) 



dt-i ' 4rj_ 

where we have used ^a*^"^ = 0, ((Iia)uj = 0, and (Jlg)t^ = 0. Applying the energy average ((...))e on both sides of Eq. 
JHOl yields 



which, using in the second term partial integration and the definition of p(e), reduces to 



^ .^5(1) ^ _^55(o) . (83) 



The vanishing of the second term on the rhs of Eq. (|82|l is the result of the spin conservation of the differential 
operator V, which in turn is ensured by the identity 

(84) 



de Td{e) 



e=0 



e=0 ^/(e) 

It is crucial to note that the rhs of Eq. H83I) is independent of the fast time because a^^^ is independent of 
t_i. Integrating Eq. (|83|l with respect to t_i, the rhs therefore gives rise to a secular term, i.e. a term which is 
proportional to As a result, 53'^^'^ can be larger than SS'^'^^ for sufficiently large times. The expansion (|76|) would 
be valid only for short times, i.e. the expansion is non-uniform. Within multiple time scale perturbation theory, 
secular terms can be avoided by an appropriate choice of the time evolutions on the various time scales. The secular 
term in Eq. (|83|l can be particularly simply removed by forcing the rhs to be zero, which gives rise to the condition 

depie)—a^°'>ito,ti) = 0. (85) 







That is, a(°)(io,ti) = a^°Hti)i which, using Eq. (O, leads to 6S^°'>{ti) = a(")(ii), i.e. the time evolution of the zeroth 
order macroscopic spin polarization (and therefore the spin decay) occurs solely on the long time scale ti. Since the 
rhs of Eq. is made to vanish, we also obtain 6S'^^\t-i,to,ti) = 5S'^^\to,ti), i.e., SS'^^^ is independent of the 

fast time variable t_i. Using the definition of SS^-^^ we furthermore conclude that a'-'^) (e, tg, ii) = a'^-* ^Oj ^i)- 
Both results we need in the analysis of the 0(r]^) equations, which is necessary to determine the time evolution on 
the remaining time scale ti. 
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3. 0{ri^) equations 

To investigate the 0{ri^) equation H71|l . we substitute (|72|l into H71|l . Averaging over the angle and energy, we find 



dt 



-1 Oto Oti \ / 



where we used Sa^^'' ~ 0, = 0, (ilg)^,^^ = 0, and {T>*a'''^^)^^^ = 0. The third term on the rhs of H86|) contains 

6a'^^^ which has to be obtained from Eq. (|81|l . the angle dependent part of the 0{rp) equation. Before we proceed 
with the analysis of Eq. H86|) let us therefore turn to Eq. H81|l . Integration of Eq. H81|l with respect to the fast time 
t-i produces a secular term, {^Ija + ^g) x a*-'^^i_i, which cannot be removed, because both ^Ija + ^g and a!"^^ are 
finite. A way to avoid the resulting non-uniformity is to demand 

^ (Sa(i) = , (87) 



dt-i 

i.e., to enforce iJa*^^) (e, w, io, ^i) = ^a^'^^e.uj.tQ.ti), which reduces Eq. (|HT|l to 

[V* - -^£2 + iClc ■ C]Sa^^^'> + (diA + (ig) X SS^°'> = 
4r^ V / 

where we used in the last term d'^^\ti) = SS'^^^ti). The condition (|87|l is reminiscent of the quasi-stationarity 
assumption usually invoked in the calculation of the D'yakonov-Perel' relaxation rates. The multiple time scale 
approach enables us to identify the time scale on which this assumption holds. 

We now return to Eq. H86I) . Since Sa^^\ SS''^\ and 55'^°^ are independent of the fast time variable the whole 
rhs of Eq. (|86|l is independent of Integration with respect to i_i thus gives rise to a secular term which has to 
be removed. We force therefore the rhs of Eq. (|86|l to vanish which can be certainly accomplished if we separately 
demand 

A<55?(i) = , (89) 
oto 

-^^SS^°'> = CtLX SS^°^ + {(Qia + fig) X <5a(i)),,^ - {R)e,u.5S^°'> . (90) 
^From the first equation we find 



/ d£p(£)— aW(e,to,ii) = 
Jo oto 



(91) 



that is a^^^ is independent of to, a result which we need below. The second equation is already a precursor of the 
Bloch equation for SS'^^^ti). Although it determines SS'^^^ti) for a given Sd^^\to,ti), it is however not yet a Bloch 
equation because, at this point of the calculation, Sa^^^ is still a function of ti and to- 

To obtain a closed Bloch equation on the time scale ti alone, we now examine the to dependence of Sa^^^ {to,ti). 
Towards that end, we consider the angle dependent part of the 0{r]^) equation ifTTjl . which reads 



dt^i ' 4ti 

+ (niA + fig) X a(i) + (niA + fig) X Sa^^^ - ((f1/A + fig) X Sa^^^)^ 

- Ra(") + (R)^a(") - ^Sa^^^ . (92) 

oto 

Except for the first term, all terms on the rhs are independent of t_i, and therefore give rise to secular terms. To 
remove the secular terms, we set the undesired terms on the rhs to zero 

— ,5a(i) = (f1,^ + fig) X Ja(i) - ((f1,^ + fig) X Sa^^^)^ 

- [R - (R)^]a(") + (niA + fig) X a(i) . (93) 
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The last two terms on the rhs of this equation are independent of tg and therefore again give rise to a secular term 
if Eq. (|93|l is integrated with respect to Iq. In general, the last two terms are finite. Thus, to avoid non-uniformity 
we demand that 5a^^^ is independent of to, that is, we enforce 5a'^^ (e, w, ioi ii) = (5a(^)(e, ^i). With this constraint, 
Eq. (|93ll could be used to determine a'-^-'(£,ii) and eventually SS^^''{ti). 

Because (5a^^^ is independent of to, Eq. (|90|l is in fact a Bloch equation on the time scale ti alone. The function 
(5a(i)(£,w,ti) satisfies Eq. ijHHl), which, through 6S'^"\ti), contains ti only as a parameter. Therefore, the function 
Sa''^\e,uj,ti) instantaneously adjusts to the function SS'^-^\ti), which, in this sense, acts like a "slave field" for 
Sa'^^^{s,uj,ti). 

To make the equations determining the decay of the macroscopic spin polarization explicit, we recall t„ = 7/"t and 
go back to the original, unsealed time variable and functions. As a result Eq. H9()(l becomes a Bloch equation for the 
macroscopic spin polarization, 

^(55(") it) = Qlx ,5^(0) it) - [Fey + rM^]SS^°^ (t) , (94) 
with initial condition (55(0) = e. The Elliott- Yafet and motional-narrowing spin relaxation tensors are given by 

Tey = (R)e,a;, (95) 

TMf,SS^°Ht) = -{(njA + dg) xSa^^'>),^^ , (96) 
respectively, and Sa^^^ is obtained from Eq. (|88|) . which for convenience we state here again 

[V* - . ^ ^^^^ ^ ^ j^(o) ^ Q 

Equations ^ ~ ^ are the main result of this section. They control the time evolution of the macroscopic spin 
polarization on the long time scale, where spin relaxation, i.e., decay, occurs. Equation 1)96(1 is an implicit definition 
of the spin relaxation tensor Tmn ■ The explicit form of Tmn can be obtained by inserting the solution of Eq. H88|l , 
which is always linear in SSit), and performing the angle and energy averages. In the next section we illustrate this 
procedure for a doped quantum well subject to a small magnetic field. 

The macroscopic spin relaxation tensor contains the EY process and the motional-narrowing (DP and VG) processes. 
Due to the different angle dependences of the two main motional-narrowing spin relaxation processes, DP and VG 
processes, Eq. (|96|l splits for isotropic semiconductors into two separate terms, Fdp and Fvg- Accordingly, for 
isotropic semiconductors, a Matthiessen-type rule holds for the total spin relaxation tensor, T — Fey + Fdp + Fyn, 
and, as a consequence, for the spin relaxation rates, which are the diagonal elements of the relaxation tensors. pLOl 
The quenching of the motional-narrowing processes due to the orbital motion is contained in Eq. I|88|) through the 
term proportional to Qc- [l^^Ol The on- and off-shell relaxation rates appearing in Eq. ((88(1 are total relaxation rates 
due to whatever scattering processes are included. A Matthiessen rule holds separately for the on- and off-shell rates. 
The energy average ((...))e is defined in Eq. ((751) . Most importantly, it contains a weight function p{e) defined as the 
solution of Eq. (|551) . This function describes the energy dependence of the quasi-stationary spin polarization which 
appears on the short time scale because of fast, spin-conserving inelastic scattering processes, i.e. because of energy 
relaxation and diffusion. Once the quasi-stationary spin polarization is established, it slowly decays on the long time 
scale set by the spin non-conserving terms in the Fokker-Planck equation. 



IV. APPLICATION TO QUANTUM WELL STRUCTURES 

In this section we apply the formalism to a quantum well at low enough temperatures, where electron- impurity 
and electron-electron scattering dominate. We arc here particularly interested in the effects of Pauli blocking and 
inelasticity. For illustration, we focus therefore only on the DP process, which, for small to moderate magnetic fields, 
is usually the dominant spin relaxation process. Moreover, if the magnetic field is small enough, the time fc it takes to 
complete a cyclotron orbit is much longer than any of the intrinsic scattering times fi (i — f, d, 1.), and the quenching 
effect of the magnetic field can be ignored. 

We consider a symmetric GaAs quantum well, grown in the [001] direction, which is also the quantization axis 
for the electron spin. Due to the assumed structural symmetry, there is only bulk inversion asymmetry giving rise 
to DP spin relaxation. |40l | As in the bulk case, we treat the two states at the conduction band minimum explicitly 
and include a large set of states perturbatively, up to third order, to include the effect of bulk inversion asymmetry. 
For energies close to the band minimum, the Hamiltonian for the quantum well can be cast into the form JQ). The 
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spin ofF-diagonal term in the quantum well Hamiltonian is the bulk spin ofF-diagonal term averaged over the envelope 
function of the conduction subband. Assuming for simplicity infinite confinement and restricting the calculation to 
the lowest conduction subband, we find (neglecting cubic terms in k) 



COS( 



hnf^is,^) = 2<5o^/i(J)2 I sine 



TiA{e) 



KiA{<jy) , (97) 



where we have defined a precession rate 1/tia{,£) = Cf^^ with Cf^ = 25q{'k / . 

Since we are only interested in the DP spin relaxation tensor, we neglect in Eq. (|88|l the torque force due to 
f2g. Because we furthermore assume small magnetic fields, we also ignore the orbital motion of the electrons. The 
separation Ansatz, 

5a(i)(e,0,i) = T{e)KiA{(t>) x S^^\t) , (98) 
then reduces Eq. H88|l to a scalar differential equation, 

4T_i(e) TiAye) 

which determines the generalized relaxation time T(e). Because the differential operator V* accounts for inelastic 
scattering, we have to conclude that even on the long time scale, where the spin polarization decays, inelasticity 
cannot be ignored. Thus, inelastic scattering processes not only determine the initial condition for the decay stage 
but they directly affect the time evolution (of the macroscopic spin polarization) in the decay stage. Multiplying from 
the left Eq. H99|l by p(£)t(£), integrating the resulting equation over e, and using condition l|84|) yields an equivalent 
differential equation, 

P - -l—\p{e)T{e) + = , (100) 

which can be also used to determine T(e). 

Inserting Eq. (|98|l into (|96|l , ignoring the fig term, and performing the angle and energy averages finally yields for 
the DP spin relaxation tensor 



1 
1 
^DP \ 2 



Tdp = — I 1 I , (101) 



with the DP spin relaxation rate given by 



— ^A — )e. (102) 

tdp tia 



where the energy average is defined in Eq. H75|) . To determine the fimction p(£), we integrate Eq. H58|) . which gives 

p{e) = , (103) 



where we used again condition H84|l and introduced the dynamical friction and diffusion coefffcients in e-space, v{e) — 
e/Tf{e) and ^(e) = e^/T(i(e), respectively. Integrating once more, we obtain 

p{e) = p(0)e^« ^ "^-f , (104) 
with w' {e) = dw{e)/de and a normalization constant ^(0), which we fix according to 

/ d<j) dep{e) = 1 . (105) 
Jo Jo 
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FIG. 3: Distribution function (for the equilibrated state at t — > <x), dynamical friction coefficient, dynamical diffusion coefficient, 
and angle randomization coefficient (from top to bottom) for a modulation doped quantum well at an electron density n — 
4 X lO^cm"^ and T = lOK. For the calculation of the coefficients, only electron-electron scattering is taken into account. The 
unit of £ is 14.93 meV. The units of the functions v, w, and u, are 22.68 meV/ps, 338.44 (meV)^/ps, and 1.52 1/ps, respectively. 



Note, because of condition (|84|l . the integral in Eq. H1U4|I is zero for e ^ 0. That is, p{e) is well-defined for e ^ 0. 

The Bloch equation (|94|1 for the macroscopic spin polarization has to be solved with the spin relaxation tensor H1U1|) 
and taking the particular geometry of the experimental set-up into account. Here, we consider the case of Kerr or 
Faraday rotation experiments, where the small magnetic field, which causes the spin precession, is along the x-axis. 
The propagation direction of the pump and probe pulses is assumed to be perpendicular to the quantum well plane, 
i.e., parallel to the z-axis (growth axis). The initial spin polarization is therefore along the z-axis, i.e. e = (0, 0, 1)^, 
and the probe pulse monitors the decay of a spin polarization which precesses in the yz-plane. Note that the spin 
decay in the yz-plane plane is not isotropic (Fyy ^ T^z). Assuming tdp ^ the solution of the Bloch equation is 

SS{t) = -sinnLt e"", (106) 

y cos flLt J 

where the decay rate of the spin polarization is given by the arithmetic mean of the decay rates in y- and in z-direction: 

r = (ryy + r,,)/2-3/(2TDp). 

The results presented below are for a L = 25nm GaAs quantum well. The parameter needed to specify h(lf^{k) 
is ^0 — O-OGfi^/ \J (2m*)^eg. (1^] The remaining parameters, such as the effective CB electron mass or the static 
dielectric constant ct, (needed for the Coulomb matrix element) can be found in standard data bases. [G^I Numerically, 
we first calculate v{e) — e/Tf{£), w{£) = £^/Td{e), which define the differential operator V, and the on-shell rate 
u{£) = l/4Tj_(e) taking electron-electron and electron- impurity scattering into account (see Appendix). We then 
determine p{e) from Eqs. Hl()4(l and H1Q5|I . Finally, we solve the differential equation (|1()0|I for T(e) numerically and 
obtain the DP relaxation rate 1 /tdp by numerically integrating Eq. Hl()2(l . 

In Figs. 121 and 0] we show, for T — lOK, the dimensionless functions v{e) — e/Tf{e), w{e) — e'^/Td{e), and 
u(e) = l/4r_L(£) for a modulation doped quantum well with electron density n = 4 x lO^cm"^ (non-degenerate 
electrons) and n = 4 x lO^^cm"^ (degenerate electrons), respectively. In a modulation doped quantum well electron- 
impurity scattering is negligible because of the spatial separation between the dopants and the electrons, we take 
therefore only electron-electron scattering into account. The physical content of the functions v{e) and w{e) is that of 
dynamical friction and diffusion coefficients (in e-space) for the "test" spin polarization resulting from the scattering 
between the spin-polarized "test" electrons and the equilibrated "field" electrons. The function u{e) denotes the 
on-shell scattering rate arising from the "test" electron's elastic scattering off "field" electrons. It randomizes the 
angle 0; hence u{s) can be interpreted as an angle randomization coefficient. 

In Figs. 13 and 0] we also show the electron distribution functions which characterize the equilibrated state at 
i — > oo. Because of our assumption that only a small portion of the total number of electrons initially contributed 
to the "test" spin polarization (i.e. 5f(k,t) <C fik)), the equilibrium distribution functions are used to determine 
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FIG. 4: Distribution function (for the equilibrated state at — > cxj), dynamical friction coefficient, dynamical diffusion coefficient, 
and angle randomization coefficient (from top to bottom) for a modulation doped quantum well at an electron density n — 
4 X lO^'^cm"^ and T = lOK. As in Fig. |3 for the calculation of the coefficients, only electron-electron scattering is considered. 
The units of e, v, w, and u are the same as in Fig. |3 

the friction, diffusion, and angle randomization coefficients (see section III.B). In other words, we approximated at 
t < oo the distributions of the friction and diffusion causing spin-balanced "field" electrons by the distributions of 
the equiUbrated electrons at t ^ oo, despite the fact that at i < c» of course a small number of electrons is still 
out of equilibrium and contributes to the finite spin polarization. Obviously, this is only permissible close enough to 
equilibrium. 

Inelastic scattering encoded in v{e) — e/Tf{e) and 'w{e) — /Td{s) gives rise to spin-conserving energy relaxation 
and diffusion, which in turn yields a quasi-stationary spin polarization (5S'st(£, </>) = p(e)e, towards which any initial 
spin polarization relaxes very quickly without losing spins. The time scale on which this spin-conserving relaxation 
takes place is set by the relaxation times Ti(£), i = /, d, _L. In Figure |S1 we show p(e) for T — lOK and three densities: 
7T, = 4 X lO^cm"^, n = 4 X IQ^^cm"^, and n = 4 x IQ^^cm"^. In the inset we again depict the corresponding 
distribution functions for the electrons which characterize the equilibrated state at t ^ oo. At very low densities, 
where the electrons are non-degenerate, p(e) is centered around £ = 0, while at high densities, where the electrons are 
degenerate, p{£) is centered around the Fermi energy for the electrons. Note, however, that p{e) describes the spin 
polarization at i = while f[e) is the equilibrium distribution of the electrons at i — > oo. 

At this point, a brief discussion about the applicability of the diffusion approximation for the linearized collision 
integral jie'' is in order. The diffusion approximation is expected to be applicable because Eq. (|37|l has the form of 
a master equation and can thus be formally expanded with respect to the momentum transfer q. For a sufficiently 
rapidly decaying transition probability W^'^{k;q), the expansion can then be truncated after the second order term 
as in Eq. (|45|l . The validity of the diffusion approximation depends therefore on the transition probability, which in 
dimensionless form reads [see Eq. (|A5ll in the Appendix] 

/oo 2 
du;\M{q)\^lfik'-u;)+n{-u;)]lmx{q,Cj)S{k-q^^-'^) , (107) 

with M[q) — l/(fZ + qs) the Coulomb matrix element, q^ the Thomas-Fermi screening wave number, and f{x) and 
n{x) the Fermi and Bose functions, respectively. Note, in Eq. H107f) we do not expand the distribution functions 
and energies. As a result, cut-off problems are avoided and the friction and diffusion coefficients, A1'^{k) and Bfj{k) 
respectively, can be calculated without restricting the (f-integration. Thus, the diffusion approximation does not ignore 
hard scattering processes with large momentum transfer; it only treats them approximately, whereas soft scattering 
processes are treated exactly. Moreover, from Eq. p07|) it follows that jie'' is actually a "phonon-type" collision 
integral, the role of phonons being played by the collective excitations of the spin-balanced "field" electrons. Indeed, 
replacing M{q) by the electron-phonon matrix element and lm.x{q, ^j-i) by the phonon spectral function, Eq. (|107|l gives 
the transition probability for spin-polarized "test" electrons scattering off equilibrated phonons. This analogy already 
suggests that a diffusion approximation is applicable to ji? . To demonstrate its validity, it is however necessary to 
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FIG. 5: Weight function for a modulation doped quantum well at T = lOK and three electron densities n — 4x lO^cm"'^ (solid 
line), n = 4 X lO^^cm"^ (dotted line), n = lO^'^cm"'^ (dashed line), taking only electron-electron scattering into account. The 
inset shows the corresponding electron distribution functions in the equilibrated state at f — * oo. 




FIG. 6: The upper panel shows P(fctyp; q) for T=10K and two electron densities: n = 4 x lO^cm"^ (non-degenerate electrons) 
and n = 4 X lO^^cm"'^ (degenerate electrons). For non-degenerate electrons, the typical momentum fctyp corresponds to the 
thermal energy (vertical solid line), whereas for degenerate electrons, fctyp is the Fermi wave number (vertical dashed line). In 
the lower panel we present the cumulant C(fctyp; g), defined in Eq. I1L)9II . which is a measure of the relative importance of the 
scattering processes with a momentum transfer less than q. 



show that for typical values of k the transition probability W^^{k; q) indeed decays sufficiently rapidly for large q. 
For that purpose we investigate 

P(fctyp; <7) = Jdcl^q W'^ihyp-, q) , (108) 

which is essentially the angle-averaged zeroth order moment of W'^'^{k — fctyp; q) with (j) the angle between fctyp and q. 
[The first and second moments appear in the calculation of the friction and diffusion coefficients Af{k) and Bfj(k)] 
For degenerate electrons fctyp is the Fermi wave number kp while for non-degenerate electrons fctyp is the wave number 
corresponding to the thermal energy. In Fig. Elwe show for the two parameter sets used in Figs. Oland^ respectively. 
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FIG. 7: D'yakonov-Perel' spin lifetime due to electron-electron scattering for a 25nm modulation doped quantum well as a 
function of electron density at three temperatures T = lOK, T = 20K, and T = 30K. The solid dots are experimental data at 
T = lOK from Ref. Isi. 



the zeroth order moment P(ktyp;q) together with the cumulant 

Jo drP{ktyp;r) 

from which we can estimate the relative importance of scattering processes with momentum transfer up to q. In the 
upper panel of Fig. we see that P{ktyp;q) indeed decays with increasing q. In the degenerate regime q — 2kp 
scattering is clearly visible, but one order of magnitude less probable than the g — > scattering process. From the 
cumulants, displayed in the lower panel, we infer moreover that in the non-degenerate as well as degenerate regime, soft 
scattering processes with q < fctyp give the main contribution (around 80-90 %). In both cases we expect therefore the 
diffusion approximation to produce reasonable quantitative results for the DP spin relaxation time. To estimate the 
error precisely is complicated. It would require a detailed investigation of the full momentum dependence of W'^'^{k] q) 
together with a reference calculation which does not invoke the diffusion approximation. The spin relaxation times 
we obtain as a function of temperature and density compare favorably with experimental results, indicating that the 
modeling of scattering processes within a diffusion approximation is sufficient for the calculation of DP spin relaxation 
times. 

We now turn to the numerical results for the DP spin lifetime tj^p . Figure {7\ shows top as a function of electron 
density for a modulation doped quantum well at three temperatures T = lOK, 20K, and 30K. Since in a modulation 
doped quantum well electron-impurity scattering is negligible, we take only electron-electron scattering into account. 
For fixed temperature, the spin lifetime first increases with electron density, reaches a maximum, and then decreases 
again. The non-monotonic density dependence of the spin lifetime follows the density dependence of the electron- 
electron scattering rate. At low densities, the scattering rate is small because of lack of scattering partners, while 
at high densities the scattering rate is suppressed because of efficient Pauli blocking. At intermediate densities, 
where the cross-over from non-degenerate to degenerate electrons occurs, the electron-electron scattering rate, and 
thus the DP spin lifetime, is maximal. The position of the maximum shifts with decreasing temperature to lower 
densities because the density, where the cross-over from a non-degenerate to a degenerate electron gas takes place, 
decreases with temperature. The relaxation time of photo-currents in optically pumped semiconductors shows a 
similar non-monotonic density dependence. |65| 

For a fixed electron density, the DP spin lifetime decreases with temperature in the low density regime and increases 
with temperature in the high density regime. The latter is because of the temperature induced reduction of the Pauli 
blocking, giving rise to an increasing electron-electron scattering rate and therefore to an increasing DP spin lifetime. 
In the low density regime, on the other hand, increasing temperature broadens the electron distribution function, 
i.e. electrons occupy states higher up in the band. The average thermal energy increases therefore and the spin 
decay occurs preferentially from states higher up in the band, where the torque force induced by the bulk inversion 
asymmetry is larger. As a consequence, the DP spin lifetime decreases with temperature in the low density regime. 



23 



1.00 



0.10 - 



0.01 




100 



electron density [10 cm ] 



FIG. 8: D'yakonov-Perel' spin lifetime for a 25nm quantum well at T = 40K and three values of ionized (donor and acceptor) 
impurity concentrations: no impurities (x — 0), ionized impurity concentration equal to the electron density {x — 1), and 
ionized impurity concentration equal to four times the electron density (a; = 4). For a; = 0, only electron-electron scattering 
contributes to the spin lifetime, whereas for x both electron-electron and electron-impurity scattering determine the spin 
lifetime. 



In Fig. [7|we also plot experimental data for T = lOK from Ref. [SC 
the agreement between theory and experiment is quite good, given the fact that our calculation is based on an 
idealized quantum well with infinitely high confinement potential. In this density regime, we expect our results 
to even underestimate the spin lifetimes, because the model for the electronic structure of the quantum well most 
probably overestimates hilf^. Indeed, the constant defining the magnitude of the splitting of the conduction subband 
(7r/L)2 ~ E"^, with E^^ the confinement energy of the lowest conduction subband. For a finite confinement 



QW 

lA 



potential, Ei is smaller than E^ — / 2m* L"^ , giving rise to a smaller splitting and, consequently, to a larger tdp- 
To obtain in this density regime better agreement between experimentally measured and theoretically calculated spin 
lifetimes an improved electronic structure calculation is clearly necessary. At lower densities, on the other hand, 
electrons are most likely localized to donors (at T — lOK thermal ionization is negligibly small) and our theory, which 
is based on a band picture, does not apply. 

In quantum wells which are not modulation doped, electron-impurity scattering due to donors and acceptors (in 
compensated samples) provides an additional, very efficient scattering process. The DP spin lifetime increases with 
scattering rate. As a result, we expect the spin lifetimes in quantum wells that are not modulation doped to be 
substantially longer than in modulation doped quantum wells. This can be seen in Fig. |S1 where we plot the electron 
density dependence of the DP spin lifetime at T = 40K for x = Q (modulation doped), x = 1 (uncompensated 
quantum well with equal impurity and electron density) , and x = 4 (compensated quantum well with impurity (donor 
and acceptor) density four times the electron density). As expected, the spin lifetimes increase with x for all electron 
densities. The increase is however not uniform, with the largest increase taking place at high electron densities, 
where Pauli blocking very effectively suppressed the DP spin lifetime in the modulation doped quantum well. The 
electron-impurity scattering rate is not affected by Pauli blocking and leads therefore to a substantial enhancement 
of the DP spin lifetime at high doping levels. 

For a fixed density the character of the electron gas also changes with temperature. In particular, increasing 
temperature pushes the electron gas into the non-degenerate regime. For some temperature, the cross-over from a 
degenerate to a non-degenerate electron gas occurs, electron-electron scattering is particularly strong, and we expect 
spin lifetimes to be enhanced. 48] This effect is demonstrated in Fig. |51 where we plot for an electron density 
n = 3 X lO^'^cm"^ the temperature dependence of the spin lifetime for a; = 0, 1, and 4. In the modulation doped case 
{x — 0) the enhancement of the spin lifetime in the temperature range where the cross-over from degenerate to non- 
degenerate takes place can be most clearly seen. For finite x, the spin lifetimes are for all temperatures longer than 
for X = Q. The enhancement is again not uniform. At high temperatures it is very small. While at low temperatures 
it is very large, because in that range the Pauli blocking leads to a strong suppression of the DP spin relaxation rate 
in modulation doped quantum wells. In fact, the spin lifetimes at low temperatures saturate by a value set by the 
electron-impurity scattering rate. The maximum in the spin lifetime is therefore less pronounced (or even disappears 
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FIG. 9: D'yakonov-Perel' spin lifetime for a 25nm quantum well as a function of temperature at an electron density n = 
3 X lO^^cm"^ and three values of ionized (donor and acceptor) impurity concentrations: no impurities (a; = 0), ionized impurity 
concentration equal to the electron density (a; = 1), and ionized impurity concentration equal to four times the electron density 
(a; = 4). As in Fig. |S| for x — 0, only electron-electron scattering contributes to the spin lifetime, whereas for x ^ both 
electron-electron and electron-impurity scattering determine the spin lifetime. 
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FIG. 10: D'yakonov-Perel' spin lifetime for a lOnm modulation doped GaAs quantum well as a function of temperature at an 
electron density n = 1.86 x lO^'^cm"^ (solid line). The solid boxes are experimental data from Ref. |4^. 



completely) in quantum wells that are not modulation doped. 

In Fig. IIUI we finally compare for a lOnm modulation doped GaAs quantum well with an electron density 
n — 1.86 X lO^^cm"^ the theoretically obtained temperature dependence of the D'yakonov-Perel' spin lifetime with 
the experimentally measured temperature dependence of the D'yakonov-Perel' spin lifetime. Below 50K the agree- 
ment between the experimental data points and the theoretical results is quite reasonable, suggesting that in this 
temperature range electron-electron scattering is the main source of electron spin relaxation in these samples. This 
conclusion is also supported by the theoretical results obtained by Glazov and coworkers. 48] In contrast to the the- 
oretically predicted non-monotonic behavior with a maximum at T w 50K, the experimental results suggest that, for 
this electron density, the D'yakonov-Perel' spin lifetime t^p grows monotonically with temperature. This is probably 
due to the neglect of phonons in the theoretical modeling. The maximum of tdp is a consequence of electron-electron 
scattering. Additional scattering processes destroy the maximum. For instance, electron-impurity scattering increases 
the spin lifetimes at low temperatures, resulting in a monotonically decreasing spin lifetime (see a; ^ curves in Fig. 
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The samples in Ref. Wi] are high-quality quantum wells, where electron-impurity scattering should be negligible. 
The electron density n = 1.86 x lO^'^^cm^^ is however rather high. As a result, the temperature, at which tjjp is 
expected to be maximal, falls in a temperature range, where electron-phonon scattering is significant. We expect a 
calculation, which takes electron-electron and electron-phonon scattering into account, to produce a monotonically 
increasing Tjjp. 

In this section we illustrated our semiclassical kinetic theory of spin relaxation by calculating the DP spin lifetime 
for an (idealized) quantum well, at temperatures and electron densities, where electron-electron and electron-impurity 
scattering dominate. Electron-electron scattering has been chosen to illustrate the effects of inelasticity and Pauli 
blocking. Electron spin lifetimes due to electron-electron scattering have been also calculated in Refs. |3a. |44L lio. lia. 
|6^ . using however different approaches and focusing mostly on different aspects. In particular, the non-monotonous 
temperature and density dependence has not been addressed until quite recently. |4^|48] In modulation doped quantum 
wells, spin lifetimes turn out to be particularly long for electron densities and temperatures, where the cross-over from 
the non-degenerate to degenerate regime occurs. In this regime, many-body effects beyond the Born approximation 
are most probably important and should be included in a more quantitative calculation of the spin relaxation times. 
We expect however our main conclusions to be independent of the particular modeling of the Coulomb interaction. 



V. CONCLUSIONS 



Starting from the full quantum kinetic equations for the electron Green functions we derived a (semiclassical) Fokker- 
Planck equation for the non-equilibrium spin polarization, assuming small spin polarizations and soft scattering. The 
Fokker-Planck equation conceptualizes the non-equilibrium spin dynamics in terms of a "test" spin polarization, 
comprising a small number of spin polarized "test" electrons, which scatter off an equilibrated bath consisting of 
impurites, phonons, and spin-balanced "field" electrons. Because of the scattering, the bath causes for the spin 
polarization dynamical friction, diffusion, and relaxation (decay) . We then empolyed a multiple time scale perturbation 
approach to separate the fast spin-conserving from the slow spin non-conserving time evolution. As a result, we 
extracted from the Fokker-Planck equation a Bloch equation which controls the time evolution of the macroscopic 
(fc-averaged) spin polarization on the long time scale, where the spin polarization decays. Our semiclassical approach 
accounts for elastic and inelastic scattering and avoids the ad hoc energy averaging of on-shell spin relaxation rates. 
Instead we show that the weight function is intimately linked to the "quasi-stationary" spin polarization, which is 
the terminating state of the fast, spin-conserving time evolution taking place immediately after spin injection. The 
diagonal elements of the macroscopic (/c-averaged) spin relaxation tensor are the spin lifetimes. They are either given 
by an energy averaged spin-flip rate (EY process) or an energy average of a generalized relaxation time multiplied by 
a precession rate (DP and VG processes). 

The formal development of our approach is based on a generic model for non-magnetic III-V semiconductors and 
treats EY and motional narrowing (DP and VG) spin relaxation processes on an equal footing. We also allowed for 
orbital motion of the electrons in a strong magnetic field, which potentially leads to a quenching of the motional 
narrowing-type spin relaxation processes. The derivation of the Fokker-Planck equation is independent of dimen- 
sionality and, as long as a soft scattering regime can be identified, also of the scattering processes, which enter the 
Fokker-Planck equation in the form of dynamical friction and diffusion coefffcients, which have to be worked out 
separately for each scattering process. 

To illustrate our formalism we applied it to a quantum well at low temperatures, where electron-electron and 
electron-impurity scattering dominate. We explicitly constructed the friction, diffusion and angle randomization 
coefficients characterizing the (symmetry-adapted) Fokker-Planck equation for that particular situation and calculated 
the DP spin lifetime at vanishingly small magnetic field as a function of electron density and temperature. We found 
that for fixed temperature (density) the density (temperature) dependence is non-monotonic. Spin lifetimes are 
particularly long for densities and temperatures, where the cross-over from a non-degenerate to a degenerate electron 
gas occurs. Spin lifetimes in compensated quantum wells are always longer than in modulation doped quantum 
wells with the same electron density. The enhancement of the spin lifetime is particularly strong for densities and 
temperatures where Pauli blocking is most efhcient in suppressing the DP spin lifetime (due to electron-electron 
scattering) in modulation doped quantum wells. 

Various extensions of our approach are conceivable and constitute research directions for the future. Semiconductor 
structures with structural inversion asymmetry [s^l and/or native interface asymmetry 31] can be studied within 
our approach by augmenting the model Hamiltonian by the corresponding spin off-diagonal Hamiltonian matrix 
elements. In particular, the role of the linear collision integral j'^\f,5S], which does not affect spin lifetimes in 
isotropic semiconductors, should be reinvestigated, e.g., for an asymmetric quantum well where spin lifetimes can be 
particularly long because motional narrowing processes due to bulk and structural inversion asymmetry can be made to 
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cancel each other. '43*, "4^ Since /g\f,SS] potentially mixes spin relaxation channels, it could affect the cancellation. 
A Fokker-Planck equation of the form (|27|1 . perhaps augmented by additional driving terms, could be the starting 
point for a systematic calculation of spin transport coefficients (e.g., spin diffusion length) for spatially inhomogeneous 
systems, such as interfaces or biased heterostructures. Finally, non-linear effects due to large spin polarizations could 
be studied either at the level of the matrix-Boltzmann equation for the electronic density matrix ^51. ,.52J or, if the 
diffusion approximation is used to simplify the collision terms, at the level of a "Fokker-Planck-Landau equation" 
for the spin polarization, where the differential operator describing spin-conserving scattering events as well as the 
spin-flip tensor explicitly depend on the spin polarization and the distribution of the spin polarized electrons. 
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APPENDIX A: CALCULATION OF RELAXATION RATES 

In this appendix we calculate for a quantum well the dynamical friction and diffusion coefficients Ai{k) and Bij(k) 
taking electron-electron and electron-impurity scattering into account. As a result, we obtain the relaxation rates 
1/Tf{e), 1/Td{e), and 1/T±{e) which define the differential operator ^{k) as well as the spin- flip tensor R(fc). Since 
we discuss in this paper only the DP process quantitatively, which originates from the interplay of the momentum 
scattering encoded in the differential operator and the torque force due to inversion asymmetry, it suffices to 

give explicit expressions only for the differential operator 2?(fc). The derivation of the spin-flip tensor R(A:) proceeds 
along the same lines. 

Within the diffusion approximation, the spin conserving collision integral jj)^'' [/, 6S] becomes a Fokker-Planck 
differential operator H45|l with dynamical friction and diffusion coefficients, A'^ and B^j, defined in Eqs. (|46|l and H47|l. 
respectively. First, we consider electron-electron scattering and calculate Af^ and Bfj. To avoid the cut-off problem 
at large momentum transfers, which usually plagues the diffusion approximation to electron-electron scattering, we 
keep the full integrands in Eqs. 146|) and (|47|l . i.e., we do not expand the distribution functions and energies with 
respect to the momentum transfer q. Using the identities f{k)f{k± q) — [f{k ± q) — f{k)]n{e{k) — e(fc ± q)) and 
(5(e(fc) -I- e(fc') — e{k — q) — e{k' + q)) ~ J dLLj5{e{k) — e{k — q) ~ lo) ■ 6{£{k') — e{k' + q) + oj), we rewrite the transition 
probability l|lT|l into 

/oo 
du^\V{q)\Hv^x{q.^)H-^) + f(k-q)]5{e{k)-e{k-q)-uj) , (Al) 

where we introduced the susceptibiHty of nomntcracting electrons 

Imx(<z,^) = 27r^[/(fc + q)-/(fc)]^(e(fc + q)-£(fc)-a;) . (A2) 

k 

The Coulomb potential V{q) is taken to be statically screened with the screening length given by the Thomas-Fermi 
expression. Had we allowed for dynamical screening, V{ci) V{q,uj) = Vo{q)/e{q,uj). If e{q,oj) is approximated by 
the RPA expression, the resulting Fokker-Planck equation would be at the level of a quantum analog to the Lenard- 
Balescu equation. |56i | The calculation of the relaxation rates presented below could be also performed with this more 
general expression for the Coulomb matrix element. For simplicity we present here however only the results for the 
statically screened Coulomb potential. 

To proceed, we introduce dimensionless quantities, measuring energies and lengths in scaled Rydbergs and Bohr 
radii, respectively. In particular, we use Rq = Rq/s, oq = Vsaoj with RqUq = h?/2mo and — 2y/sRodo, and choose 
s such that Rq = ImeV. The dimensionless Fokker-Planck operator has the same form as in Eq. (|45|l with dynamical 
friction and diffusion coefficients given by 

Ar{k) = C^e J dqq.,W{k-q) , (A3) 

Bt!ik) ^ ^ j dqq,q,W{k-q) , (A4) 
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with = sm* / el-Kirifj and 

J-oo ('7 + 9s)^ 2 2 

where is the Thomas- Fermi screening wave number and £b is the static dielectric constant. The function F(k'^, uj, q) 
originates from the statistics of the electron gas and is given by 

F{e,uj,q) = lmx{q,Lo)N{k,Lo) , (A6) 

/2 
dfc[/(fc2 +co)~ fi^mk . g + |- - |) , (A7) 

N{k,uj) = /(fc2 - w) , (A8) 

with f{x) and n{x) the Fermi and Bose functions, respectively. To calculate Af'^(k) and Bf^{k) for a [001] quantum 

well, we first evaluate the integrals for a fixed coordinate system in which k = kcx and then rotate to an arbitrary 
coordinate system. The result can be cast into the form 

A-(fc) = -^G(fc)fc, , (A9) 
2 2 

respectively, with three functions defined by 



Bf;(fc) = -H{k)S,, + —E{k)hkj , (AID) 



= -5 /'"^"^ dq -^f^^'-f) - , (All) 

2 i-oo (9 + 9.)' 



H(fc) d^i dq ' VT^^ , (A12) 



""^'^ ^ --j^j^L. ' VI ■ ^^''^ 

Note that these integrals are well defined. The range of integration originates from the (/)-integration which also gives 
rise to the factor involving z — {q"^ + uj)/2kq. With Eqs. (|A9|) and (|A10|I and a transformation to the radial variable 
e — k"^ , the spin-conserving part of the dimensionless electron-electron collision integral is given by (recall that in two 
dimensions SS{e,4',t) contains the factor J{e)/ {2Tr)'^ns = I/Stt'^Us) 

Ji°J[f,SS]^V,e{k)SS{e,<l>,t) , (A14) 

with 

T>ee{k) = -— Wee(e) + ^^ee(e) " Uee{e)C'^ . (A15) 

Here we have introduced the total angular momentum operator in two dimensions, C = —id/dcf), and the friction, 
diffusion and angle randomization coefficients, Ved^), Weei^), and Ueds), which are linear combinations of the fimctions 
G(fc), H{k), and E{k) taken at fc = ^/e. Specifically, they read: 

2C™ r ^ z-?"- ujF{e,u,q) 1 
= ^Gll^e d,^ r^^dq f ^ , (A17) 



— oo 



uUe) = ^ dco dq^P^^^^/l^ . (A18) 



the sign of uj, qmax = Q^axi^ 

(e,w) = -Ve + Ve^, < , (A19) 



The range of integration depends on the sign of uj, qmax = Qmaxi^^ + '?maa;(^i ^)®('^) 9' 

9™«(£'^^)0(-^) + 9m™('^'^)0(^)' with 





— UJ 


'e 


— UJ 


'e 


— UJ 



q^axie,^) = V£ + y e-uj , ui < , (A20) 
^) - - , < < e , (A21) 



Qmaxis,^) ^ Vi + V^^, 0<Lo<e. (A22) 
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The functions Vee{£) and Wee{£) are the dynamical friction and diffusion coefficients for the spin polarization in e- 
space. They originate from the scattering of the "test" electrons comprising the spin polarization with the equilibrated, 
spin-balanced "field" electrons. Because the scattering is inelastic, the "test" electrons gain {lo < 0) or lose {uj > 0) 
energy by scattering off "field" electrons. The on-shell function Mee(e) describes randomization of the angle ip. 
The integrals defining t;ee(e), Wee{s), and Ueei^) have to be done numerically. The singularities are integrable and 
Gaussian integration proved to be efficient. The limiting values are lime_»oWee(e) — uo/e, lime^oWeeis) — 0, and 
\m\^QVeeis) = vo- Moreover, Weeis) = vqE for e — > 0, i.e. dwee{0) / ds = Wee(O), which is essential to guarantee spin 
conservation of the differential operator I?ee(fc) [cp. Eq. (|84ll ] . The structure of the differential operator suggests to 
write = e/tJ'^, Wee — ^"^ / ^ "ee — ^ / 4:t'^ , with ^ / T^'^ , ^ / T^j^^ , and I/t]^ relaxation rates describing energy 
relaxation, diffusion, and randomization of the angle due electron-electron scattering, respectively. 

The calculation of the (on-shell) relaxation rate due to electron-impurity scattering proceeds along the same lines. 
The starting point is Eq. (|37|l . specialized to electron-impurity scattering, that is, with W'^{k,q) given by Eq. I|40|l . 
Going through the same steps as in the case of electron-electron scattering yields for the spin-conserving part of the 
dimensionless electron-impurity collision integral (as before SS{e,(f>,t) contains the factor l/87r^ris) 

4°\f, SS] - V,,{k)SS{e, 0, t) , (A23) 
where the differential operator is now given by 

Ve^ik) = -ue,(e)£2 , (A24) 

with 

^e^{e)^^ dq . I ,, ^l-z^ , (A25) 

C™ — s{4:TTm* niUQ / elnio) and z = q/2y/e, where rii is the sheet density of the impurities. Because of the elasticity 
of electron-impurity scattering, the differential operator contains only an on-shell term. The function Uei = I/AtJ^ 
defines the relaxation time due to electron-impurity scattering. It only causes randomization of the angle 0. The total 
spin-conserving collision integral, taking electron-electron and electron- impurity scattering into account, is given by 
D{k)=D,,ik)+D,,ik). 

Similar expressions can be derived for electron-phonon scattering. For bulk semiconductors the calculation proceeds 
along the same lines with the obvious modifications due to the additional angle integration. 
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